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Abstract: Series for the Wilson functions of an "environmentally friendly" renor- 
malization group are computed to two loops, for an 0{N) vector model, in terms of the 
"floating coupling", and resummed by the Pade method to yield crossover exponents for 
finite size and quantum systems. The resulting effective exponents obey all scaling laws, 
including hyperscaling in terms of an effective dimensionality, d^jf = 4 — 'jj^, which repre- 
sents the crossover in the leading irrelevant operator, and are in excellent agreement with 
known results. 
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Physical systems can exhibit different types of scaling behaviour in different asymptotic 
regimes. The crossover between such asymptotic regimes is important both theoretically, 
and experimentally. One may think of a crossover as being induced by some "environ- 
mental" variable. Two of the most interesting crossovers are: those induced by finite 
size effects [1], and those induced by quantum effects [2], [3]. For finite size systems the 
environmental variable is L, the system size, while in quantum systems it is the inverse 
absolute temperature, Ph. Other environmental factors that may induce a crossover are: 
long range interactions, anisotropic interactions, external fields, boundary conditions, etc.. 

The main difficulty in treating systems which exhibit a crossover is that the quali- 
tative nature of the effective degrees of freedom (DOF), i.e. the fluctuations, changes 
significantly as a function of scale, being very sensitive to the environment in the crossover 
region. The renormalization group (RG) is our most powerful tool for investigating how 
physical systems change as a function of scale. If one views the RG transformation as a 
"coarse graining" procedure, one must ask whether a particular coarse graining is captur- 
ing the qualitative changes associated with the crossover. We call an RG which tracks 
the changing nature of the effective DOF — "environmentally friendly" . Momentum shell 
integration is, in principle, a naturally environmentally friendly form of renormalization, 
however, the associated RG's are not easily computed beyond lowest order in perturbation 
theory. In addition one must be careful not to throw away possible important environ- 
ment dependence by arguing that integrating out large momenta k ^ g, where g is the 
characteristic scale set by the environment, should be g independent, as this leads to RG 
flow equations which when propagated to scales k < g coarse grain effective DOF which 
are a very poor representation of the system's fluctuations at that scale. 

The most accurate RG results for properties of critical systems, such as critical ex- 
ponents, have been achieved by applying field theoretic techniques [4], and are in very 
impressive agreement with experiment [5]. Three approaches have been used: s expan- 
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sion [6], expansion [7], and fixed dimension perturbation theory [8]. The first fails in 
crossovers where the upper critical dimension changes, such as in finite size crossover. The 
second fails in crossovers where the order parameter can change its symmetry, such as 
bicritical crossover. We adopt the spirit of the fixed dimension approach, if not the letter, 
in the context of environmentally friendly renormalization. 

Field theory historically, as distinct from momentum shell integration, has empha- 
sized the role of ultraviolet (UV) divergences, which are independent of infrared scales and 
therefore environment insensitive. The result is an RG which does not track the changing 
nature of the effective DOF, and this leads, typically, to a breakdown of perturbation the- 
ory. More environmentally friendly RGs have been implemented in some contexts. Amit 
and Goldscmidt [9] introduced the concept of generalized minimal subtraction (GMS) 
when considering crossover at a bicritical point. Their results for 7e//, however, differ from 
those found using momentum shell integration [10], the latter find a characteristic "dip" 
in the effective exponent curves. Our methodology applied to a bicritical crossover [11] 
gives results in agreement with the momentum shell approach. GMS was also applied to 
uniaxial dipolar ferromagnets in [12], the results of our analysis [13, 11] differ somewhat. 
Differences can be understood, since our renormalization procedure puts complete Feyn- 
man diagrams into the Wilson functions, GMS does not. When the RG equation is solved 
all the diagrammatic information in these functions is "exponentiated", the remainder 
must be taken into account perturbatively, in some way. Other related work is that of 
Schmeltzer [14] who calculated 7e// to one loop for a three dimensional quantum ferroelec- 
tric and Lawrie [15] who considered dimensional crossover for d-dimensional quantal and 
d + 1 dimensional finite-sized Ising models for 3 < d < 4 using an s expansion. Unlike 
our method the e expansion could not capture the crossover between two non-trivial fixed 
points as the upper critical dimension changes across the crossover. Field theoretic results 
for dimensional crossover in a fully finite geometry or a cylinder have been obtained [16] 
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but the techniques used have not been extended to the case of a system with more than 
one fixed point. 

Though our general approach is apphcable to a very wide class of crossovers [11,17], we 
restrict our attention to finite size crossover and quantum/ classical crossover. We begin 
with the "microscopic" Landau-Ginzburg- Wilson Hamiltonian 

H[iPB\ = ^ j d'^x ^(V^b)2 + \ml^l + \tB{x)^l + - Hb{x)^b (1) 

which describes either: a layered d+1 dimensional system, of thickness L; or a (i dimen- 
sional quantum system, with L = ph, P being the inverse temperature. We will assume 
the order parameter possesses an 0{N) symmetry, the case AT = 1 of quantum/classical 
crossover represents an Ising model in a transverse magnetic field. In the finite size case 
+ tB = T — Tq, and in the Ising model in a transverse field + = F — Fq. Here, 
Tq and Fg are the critical temperature and transverse field respectively, in the mean field 
approximation. 

An L dependent renormalization is necessary to obtain the desired environmentally 
friendly RG. Use of L dependent normalization conditions achieves this and ensures that 
all the Feynman diagrammatic information is exponentiated in the solution of the resulting 
RG equation. The relation between the bare and renormalized vertex functions is r^'^^ = 

_N_ 

^ Z~2^r(-^'-^) The renormalized dimensionful coupling is similarly related to the bare 
one by = Z'^^X (see [4] for a discussion of the notation). We choose the normalization 
conditions 

(i) r(2)(fc = 0,t = Ac2,A,L,Av) = Av2 (ii) f^(/c,t = «2^A,L,Av)|,=o=l 

(iii) F(4)(/c = 0,t = K2,A,L,K) = A (iv) T('^^'^\k,t = h?, \ L, k) = 1 

which specify Z^, Z^2 and Z^- Condition (i), together with the multiplicative renormal- 
ization of t, implies that t is proportional to T — Tc{L) for the finite size system, and 
F — Fc(/3) for the quantal Ising model, i.e. that one is measuring temperature/field devia- 
tions relative to the L dependent critical point. We are assuming here that the system can 
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exhibit critical behaviour for any value of L, which restricts our attention to d > 1 in the 
case of A'' > 1, but in no way restricts the generality of our approach however. We have 
applied our methods successfully at one loop to dimensional crossover in a non-linear a 
model [18] and find results in qualitative agreement with those of [3] where a momentum 
shell integration approach was used, also at the one-loop level. We believe our methods are 
more easily extended to higher orders. If a normalization condition with L = oo had been 
used, then temperature/field deviations would be measured relative to Tc(oo) or Fdoo), 
the critical temperature of the bulk system, or critical field of the P = oo quantum system 
respectively. Note that we are here using an RG which runs the renormalized temper- 
ature parameter in distinction to the Callan-Symanzik equation which runs the physical 
correlation length. 

The RG equation can be viewed as a simple consequence of the fact that the bare theory 
is independent of the arbitrary renormalization scale k at which we choose to define our 
parameters, i.e. i^^^^^^ = 0. Using the relation between the bare and renormalized vertex 
functions and expressing things in terms of the renormalized parameters the infinitisimal 
form of the RG equation then becomes 

with 

1 dZ^ 1 dZ-} p(^x) 1 dZx 

^^^z^^iii^-^ = -ir = ^^ = z-x^i[^ 

The functions 7^, 7^2 and 7;^ are the Wilson functions. They are explicitly L dependent 
due to the normalization conditions (2) and all the physics of the crossover can be gleaned 
from them. 

A suitable coupling, with respect to which perturbation theory can be performed, is 
the fioating coupling [19,11], /i, which is chosen so as to make the quadratic term in I3{h) 
have unit coefiicient. Our perturbation theory is then carried out at the level of the Wilson 
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functions in terms of h. The expressions obtained are, however, only the leading terms in 
an asymptotic expansion of the functions (3{h, z), 7,^2 (/i, z) and 7(^(/i, z). We use [2,1] Pade 
approximants to resum these asymptotic series obtaining 

(3{h,z) = -e{z)h + y ^ (5) 



and 



+ h 

VC.-)- (jV + 8)l + 6^(/l(.)-|/2W)/. " 

where the functions £, /i and /2 depend on d and z — kL but are independent of N . The 
original non Pade resummed series can be recovered by expanding 1/(1 + xh) ~ 1 — xh. 
We will take the solution of (5) as our perturbation parameter. After these equations are 
solved it is then inappropriate to do any further expansion. 

The functions e{z), fl{z) and f2{z) are the basic building blocks, their specific func- 
tional form depending on the particular crossover in question. e{z) can be thought of as 
being a measure of the "effective dimensionality" of the system. The functions /i and 
/2 for general d and the crossovers of interest here can be found in [17]. For (i = 3, the 

expressions become especially simple, we find e{z) = 1 — 2;^ln(^^m~^) 

n 

r 1 r 1 - Ml 1 .. r 1 I ^ y i 

J. , ^, ^nx,ni . . .ni,n2 

/l(^) = 2 ^ 1 , and = 4- 



n n 

with mi = (1 + mi2 = (1 + + "-2)^)^ M = mi + 777,2 + mi2. We plot e{z), 

fl{z) and f2{z) against ln{l/z) in Figure 1. 

Effective critical exponents defined as logarithmic derivatives of the associated ther- 
modynamic quantities with respect to T — Tc{L) at fixed L for the finite size crossover 
and with F — Fc(/3) for fixed /3 in the quantum problem, using the above RG [11,20] 
can be shown to obey all the usual scaling relations including hyperscaling. The usual 
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dimension is replaced by the effective dimension d^ff = A — ^\ which reflects the chang- 
ing importance of the leading irrelevant operator. As a consequence these exponents 
are related to the Wilson functions through: i^e// = 1/(2 — 7^2), ?7e// = 7^, 7e// = 
(2-7,/.)/(2-702), a,ff = (7A- 27^2)7(2-7^2), = (2 - 7^ + 70)/(4 - 27^2) and 

^eff = (6 — 7a — 7^)/ (2 — 7a + 7^)- Analogous effective exponents associated with varia- 
tions with respect to L at fixed T, and (3 at fixed F can also be defined and computed. 

We present our results in graphical form in figures 2 through 5. In all graphs the 
horizontal axis is ln(^^/L), the different curves correspond to = (polymers), A = 1 
(Ising model), A = 2 (XY-model), A = 3 (Heisenberg model) and A = 00 (spherical 
model). The curves represent both a four dimensional layered geometry of thickness L 
and a three dimensional quantum model at (3 — L. The logarithmic corrections to scaling 
at the bulk end are clearly visible, their magnitude is as expected from four dimensional 
calculations. All curves are with the boundary condition /j, = 1 at ln(^^/L) = —20, the 
value of h at the initial scale parameterizes different possible crossover curves but all curves 
asymptote to the same form. In Figure 2 we plot Veffi the correlation length exponent, 
for N — 00, Ueff = ^/{deff — 2) across the entire crossover. In Figure 3 we plot ?7e//, the 
exponent which governs the fall off in critical correlations at T = Tc(L) and F = Fc(/3) for 
finite size and quantum systems respectively. This exponent is not a monotonic function 
of A but attains a maximum for some value between A = — 2 and A = 00, where it is 
identically zero. This is the least accurate of our exponents and the peak appears to be at 
A = 1, though more accurate values for this exponent suggest it occurs at higher values, 
probably A = 3. Figure 4 shows a plot of the effective specific heat exponent a^ff which 
measures how the singular part of the free energy changes as F or T varies. The extra case 
A = —2 is added here, since, in the case of dimensional crossover it is distinguishable from 
the Gaussian model due to the fact that 7^ for the latter is zero whereas for the former it is 
non-zero, being a measure of the changing effect of the leading irrelevant operator. Across 
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the entire crossover one has cte// = 2 — i/^ffdeff. Not only does one see the change in sign 
of the specific heat exponent as a function of N but one also sees that the effective specific 
heat exponent can change sign as a function of ^l/L. This is quite pronounced for the XY 
model which starts off positive, increases then turns negative at $,i ~ 100 L. It would be 
interesting, based on the Harris criterion for the relevance or irrelevance of weak disorder, 
to see whether disorder could change from being irrelevant to relevant as a function of size, 
or temperature in the case of a quantum system. In Figure 5 we plot 7;^ = 4 — deff which 
also gives information about the effective dimensionality of the system. Notice that 7;^ 
is very robust to changes in N, varying very little across the entire range of N, [—2, 00]. 
The other effective exponents can be determined from the effective exponent laws, which 
we have verified also by direct calculation. Asymptotic values of critical exponents and 
associated quantities are tabulated below. 



Asymptotic Critical Exponents 


N 


7^ 




h 


7.// 






- 2 


0^ 


0* 


1.800 


1* 


0.5^= 


0.5^= 


- 1 


0.0200 


0.145 


1.820 


1.088 


0.550 


0.351 





0.0295 


0.277 


1.785 


1.175 


0.596 


0.211 


1 


0.0329 


0.388 


1.732 


1.257 


0.639 


0.083 


2 


0.0332 


0.479 


1.675 


1.330 


0.676 


- 0.029 


3 


0.0322 


0.552 


1.621 


1.395 


0.709 


- 0.126 


4 


0.0305 


0.611 


1.573 


1.451 


0.737 


- 0.211 


00 


0* 


1* 


1* 


2* 


1* 


- 1* 



* These values are exact. 



All these values are in very good agreement with corresponding high temperature series 
and experimental results (see [21] and references therein). We believe the entire crossover 
curves are of similar accuracy. 

In this paper, using two loop Fade resummed perturbation theory for an environmen- 
tally friendly RG, we presented effective critical exponents for dimensional crossover in a 
four dimensional layered system with periodic boundary conditions and quantum to clas- 
sical crossover in three dimensions. We paid special attention to polymers, Ising model, 
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XY-model, Heisenberg model and the spherical model. Asymptotic values for the expo- 
nents of these systems were found to be in very good agreement with known results and 
experiment. Our general formalism is applicable to a wide class of crossover problems. 
Three and higher loop calculations, we believe, are quite feasible numerically by methods 
similar to that of Nickel [22]. These should yield effective critical exponents to the same 
degree of accuracy as standard critical exponents. There is merit in pursuing such cal- 
culations as our methods provide a direct and physical connection between exponents in 
different dimensions. 
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